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A progress report on experiences with a gauge fixing method proposed in LATTICE 94 is presented. In 
this algorithm, an SU(N) operator is diagonalized at each site, followed by gauge fixing the diagonal (Cartan) 
part of the links to Coulomb gauge using the residual abelian freedom. The Cartan sector of the link field 
is separated into the physical gauge field ct~J^ responsible for producing /^f rtan , the pure gauge part, lattice 
artifacts, and zero modes. The gauge transformation to the physical gauge field 

<#> is then constructed and 

performed. Compactness of the fields entails issues related to monopoles and zero modes which are addressed. 

AZPH-TH/96-18 



1. The Method 

While gauge fixing is a central tool of lattice 
simulations, the effect of lattice artifact "Gribov 
copies" remains a delicate issue, particularly in 
chiral fermion models which often rely on gauge 
fixing essentially. 

In [0 we proposed a method of gauge fixing to 
an t'Hooft like gauge which was called Gauge fix- 
ing by Hodge decomposition. We work here on a 
spatial torus in 3-dimcnsions, thus this presenta- 
tion is for Coulomb gauge (the algorithm is done 
in parallel on each time slice). Landau gauge gen- 
eralizes similarly. The method is as follows. 

• Diagonalize some operator which trans- 
forms adjointly, 0' x — G x O x G x , at each 
site. (The operator used here is the spatial 
sum of plaquette clovers and their adjoints 
at each site.) 

• Define an Abelian a M field from the links. 



1.1. Diagonalization of O 

An iterative method can be used in which the 
operator O is hit with successive SU(2) subgroup 
gauge transformations that each minimize the 
modulus of resulting off-diagonal terms. If 
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Decompose a^(x) — oqf' + + + 
Wfj,, and solve for the physical field ajL 
which minimally produces the plaquettes 
and open (monopole) strings f^ v . 
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Add the continuum zero-mode: 
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• Construct and perform the U(l) gauge 
transformation which moves the links from 
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to ajj + w 
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For O e SU(2) or su(2), only one of each 
gauge transformation is required. For SU(3) op- 
erators, about 10 iterations are required to get 

Soff-diagonal \ a ij \ < 10 ? - 

1.2. The Residual Abelian Fields 

We must define residual U(l) fields, which rep- 
resent the part of A° in the Cartan subgroup in 
the continuum limit. For SU(2) we use: 

= tan _1 (e 3 /eo) (4) 

where = e + ie k cr k . 

For SU(3) where a link has diagonal elements 
{e^e^e^ 3 }, the angles 

al = (2/3 1 -/3 2 -/3 3 )/3, and 
c£ = (2/3 2 - /3 1 - /3 3 )/3 

are suitable Cartan fields. 

2. Hodge Decomposition of a M 

Any lattice vector field a fl (x) can be uniquely 
decomposed into: 



(5) 



where <9 M (djj is the lattice forward (backward) 
derivative 

f t\ 

• otjl is the physical part of the field and so- 
ley responsible for producing 

f/tv = dftQ.^ - d v a\f\ Conveniently, 
also naturally satisfies the Landau gauge 
condition: d^a^ = 0. 

• d^ifi is the pure gauge part. 

• H,j, is the lattice harmonic part responsi- 
ble for Dirac string loops (and is the major 
source of Gribov copies). 

• is the continuum harmonic part for a 
torus, ie. a constant. 

The decomposition is formal at this stage, but 
we remark that much of the work in solving for 
the lattice a^f* is in identifying the neccessary 
parts of Hp which must be kept. 



2.1. Solving for a { J ] 

aP is determined (up to zero-modes) by 



or 



(6) 



(7) 



in Fourier space. fc M is the lattice momentum: 
2sin(7ran A j/i M ). It is thus very easy to find aj/' 
from fp V by FFT; however we must build the 
correct in stages. 

Since we want djf* to minimally reproduce the 
plaquctte angles fffl: 

. We first set /„„ - ffl. 



3. Monopoles 



For compact gauge fields, a)i must also 
the reproduce gauge invariant monopoles which 
are the ends of open lengths of Dirac string; 
thus we need to find all monopole sites on the 
dual lattice. Since the Dirac strings connecting 
these monopoles are gauge variant, we define a 
monopole-anti-monopole pairing which minimizes 
the string length between pairs. This can be done 
by simulated annealing for instance, but since this 
pairing is not unique, this step is a source of Gri- 
bov ambiguity. We could also connect pairs in 
order of finding them (violating rotational invari- 
ance), which is however fast and unique. Then, 

• We next add to / M „ = + f™ nplc 

in othcrwords, for each plaquette pierced by a 
Dirac string (as given by our minimal monopole 
pairing), we add ±2tt to / M „. 

3.1. Zero- modes of / M „: globally wrapping 
string 

Due to compactness again, may have a 
zero-mode. This zero-mode can be viewed as a 
length of non-contractable Dirac string, ie. one 
that stretches across the torus. 

^3 J dx»dx»dx»(fft + /- nplc ) =c^0 

{p,,v _L p) indicates the occurance of c p global 
non-contractable Dirac strings in the p-direction 
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in the initial configuration which must be added 

to f/ju,. 

We are at liberty to place these strings wher- 
ever we like, either randomly for pseudo transla- 
tional invariance, or along some particular axis so 
that we know where they are. Along these strings 
we add ±2tt to the perpendicular to the path 
so that = f$* + /-"P 10 + /sj bal 

We are now ready to solve eq. @, which will give 
the minimal ajj that reproduces the plaquette 
angles, monopole, and global strings, and which 
is in Landau gauge: 



3.2. One more zero-mode 

There is one more zero mode that is undeter- 
mined by eq. H which is responsible for producing 
the correct Wilson loops. After solving for a)i 
we must add to it the correct constant 

Wm = h J J dx " dxV J dx ^ (" m ~ °^ ) ) (8) 

in order that Wilson loops are preserved mod(27r). 

This is the zero-mode of the original a M field, 
and does not contribute to / M „. 

4. The Gauge Transformation 

Because of the zero-modes, the only way to 
find the gauge transformation taking us from 
— ► cuff' + is by constructing a tree, or 
in otherwords integrating the equation 

9 n9 = ( a i f ) + w fi)~ 

5. Tests 

Figure 2 shows three extremization gauge fixed 
copies derived from the initial configuration in 
figure 1. In figure 3, the only two copies ob- 
tained by the Hodge method are displayed. The 
SU(2) starting configuration is relatively smooth 
though, generated at (3 — 44, followed by a ran- 
dom gauge transformations. 




Figure 1. Initial gauge field configuration and 
string topology, on slice t = 0. The "sum of 
clovers" of an SU(2) field was first diagonalizcd 
according to the introduction, then a random 
gauge transformation was applied. 



6. Conclusions 

• The method seems to work resonably well; 
it returns a uniquely gauge fixed configu- 
ration, up to the connectivity of monopole 
pairs. However, the Cartan sector of 
SU(N) fields at typical /3 values is fairly 
rough, and thus the monopole density is 
also relatively high. 

• Closed (contractable) loops of Dirac string 
are removed, which are a primary source of 
Gribov copies in extremization methods. 

• For high monopole densities, simulated an- 
nealing seems to give poorer connectivity 
for the large number of monopole pairs than 
extremization, ie. extremization uses less 
string to connect monopoles. 

• For rough fields, \cejj \ can occasionally be 

larger than tt. This means that the phys- 

f ft 

ical (non-compact) field a]i is "clipped" 
by compactness, resulting in a site with 
9 M a M = ±2tt. 
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Figure 2. Three of the many Gribov copies from 
gauge fixing by extremization. 




Figure 3. Hodge method: only 2 final config 
tions occur. 



